REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  coilection  of  information  is  estimated  to  average  1  hour  per  response,  inciuding  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  coiiection  of  information.  Send  comments 
regarding  this  burden  estimate  or  any  other  aspect  of  this  coilection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 
Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Ariington  VA,  22202-4302. 
Respondents  shouid  be  aware  that  notwithstanding  any  other  provision  of  iaw,  no  person  shaii  be  subject  to  any  oenaity  for  failing  to  comply  with  a  coiiection 
of  information  if  it  does  not  dispiay  a  currentiy  vaiid  OMB  controi  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 
Final  Report 


1.  REPORT  DATE  (DD-MM-YYYY) 
06-01-2016 


4.  TITLE  AND  SUBTITLE 

Final  Report:  Theoretieal  Studies  of  Small-System 
Thermodynamics  in  Energetic  Materials 


3.  DATES  COVERED  (From  -  To) 
l-May-2012  -  30-Sep-2015 


5a.  CONTRACT  NUMBER 
W911NF-12-1-0146 


5b.  GRANT  NUMBER 


6.  AUTHORS 

Donaald  L.  Thompson,  Thomas  D.  Sewell 


5e.  PROGRAM  ELEMENT  NUMBER 
611102 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRESSES 

University  of  Missouri  -  Columbia 
115  Business  Loop  70W 
Mizzou  North,  Room  501 

Columbia,  MO  65211  -0001 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS 
(ES) 

U.S.  Army  Researeh  Offiee 
P.O.Box  12211 

Researeh  Triangle  Park,  NC  27709-2211 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

61404-EG.6 


12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  Publie  Release;  Distribution  Unlimited 


13.  SUPPLEMENTARY  NOTES 

The  views,  opinions  and/or  findings  eontained  in  this  report  are  those  of  the  author(s)  and  should  not  eontrued  as  an  offieial  Department 
of  the  Army  position,  poliey  or  deeision,  unless  so  designated  by  other  doeumentation. 


14.  ABSTRACT 

This  is  a  comprehensive  theoretical  research  program  to  investigate  the  fundamental  principles  of  small-system 
thermodynamics  (a.k.a.  nano  thermodynamics).  The  proposed  work  is  motivated  by  our  desire  to  better  understand 
the  fundamental  dynamics  and  thermodynamics  of  hotspot  formation  and  propagation  in  energetic  materials.  A 
better  understanding  of  systems  that  are  nanoscale  and  smaller  that  are  either  isolated  or  embedded  in  bulk  can 
provide  an  improved  understanding  of  initiation  and  reaction  propagation  criticality  (that  is,  hotspots)  in 


15.  SUBJECT  TERMS 

hotspots,  small-system  thermodynamics,  nanoenergetics,  molecular  dynamics 


— 1, 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

15.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF  PAGES 

UU 

UU 

UU 

UU 

Donald  Thompson 


19b.  TELEPHONE  NUMBER 
573-882-0051 


Standard  Form  298  (Rev  8/98) 
Prescribed  by  ANSI  Std.  Z39. 18 


Report  Title 

Final  Report:  Theoretical  Studies  of  Small-System  Thermodynamics  in  Energetic  Materials 

ABSTRACT 

This  is  a  comprehensive  theoretical  research  program  to  investigate  the  fundamental  principles  of  small-system  thermodynamics  (a.k.a. 
nano  thermodynamics).  The  proposed  work  is  motivated  by  our  desire  to  better  understand  the  fundamental  dynamics  and  thermodynamics 
of  hotspot  formation  and  propagation  in  energetic  materials.  A  better  understanding  of  systems  that  are  nanoscale  and  smaller  that  are  either 
isolated  or  embedded  in  bulk  can  provide  an  improved  understanding  of  initiation  and  reaction  propagation  criticality  (that  is,  hotspots)  in 
miniaturized  and  precisely  engineered  high-energy -density  materials.  The  proposed  work  will  improve  our  fundamental  knowledge  of 
energy  localization  and  energy  transfer  for  various  physical  situations  germane  to  traditional  and  nanoscale  energetic  formulations.  It  is 
based  on  atomic-scale  simulation  methods  -  primarily  molecular  dynamics  with  realistic  force  fields  -  and  consists  of  fundamental  studies  of 
atomic  and  molecular  clusters,  nanoparticles,  and  single-  and  multi-component  condensed  phases  for  materials  ranging  in  complexity  from 
simple  atomic  systems  through  molecular  explosives  and  interfaces  between  explosive  constituent  materials.  Exploring  how  the  fundamental 
thermodynamic  properties  and  energy  redistribution  processes  depend  on  system  size,  thermodynamic  phase,  or  system  dimensionality  will 
contribute  to  understanding  the  detailed  interactions  that  might  be  useful  in  the  design  of  advanced  energetic  formulations.  The  energetic 
compounds  of  prime  focus  will  be  nitromethane,  PETN,  TATB,  and  cyclic  nitramines;  however,  we  may  also  study  other  systems,  for 
example  inorganic  salts,  as  indicated  by  our  results  to  further  explore  fundamental  issues. 
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This  project  focused  on  investigations  of  fundamental  small-system 
thermodynamics  (aka  nano-thermodynamics)  in  organic  energetic  materials.  The 
practical  motivation  was  to  better  understand  the  fundamental  dynamics  and 
thermodynamics  of  hotspot  formation  and  propagation  of  chemistry  in  energetic 
materials.  Thus,  we  focused  on  the  behavior  of  regions  within  a  bulk  material  that 
initially  is  not  in  thermodynamic  equilibrium  with  the  surrounding  material.  We 
used  classical  molecular  dynamics  simulations  to  investigate  fundamental 
properties  at  the  atomic  level.  Two  basic  types  of  "nano-regions"  were  considered: 
Small  numbers  (including  one)  of  molecules  in  a  solid  and  interfaces  between  a  solid 
and  a  liquid. 

We  completed  three  studies;  two  of  which  have  been  published  and  the  third  in  a 
manuscript  to  be  submitted  for  publication.  The  first  study  investigated  the 
relaxation  of  an  excited  molecule  in  solid  nitromethane.  In  the  second  study 
analytical  expressions  for  the  equations  of  state  of  simple  systems  were  developed 
and  compared  to  molecular  dynamics  results.  In  the  third  study  we  characterized 
equilibrium  interfaces,  studied  the  effects  on  heat  flow  (Kapitza  resistance)  through 
interfaces,  and  the  effects  of  interfaces  on  shock  fronts.  Abstracts  of  the  articles  and 
the  manuscript  follow. 

A  molecular  dynamics  study  of  the  relaxation  of  an  excited  molecule  in  crystalline 
nitromethane,  Luis  A.  Rivera-Rivera,  Ali  Siavosh-Haghighi,  Thomas  D.  Sewell,  Donald 
L.  Thompson;  Chemical  Physics  Letters,  Vol.  608,  21  July  2014,  Pages  120-125; 
doi:10.1016/j.cplett.2014.05.065. 

Abstract:  Classical  molecular  dynamics  simulations  were  used  to  study  the  relaxation  of 
an  excited  nitromethane  molecule  in  perfect  crystalline  nitromethane  at  250  K  and  1  atm 
pressure.  The  molecule  was  instantaneously  excited  by  statistically  distributing  energy  E* 
between  25.0  kcal/mol  and  125.0  kcal/mol  among  the  21  degrees  of  freedom  of  the 
molecule.  The  relaxation  occurs  exponentially  with  time  constants  between  11.58  ps  and 
13.57  ps.  Energy  transfer  from  the  excited  molecule  to  surrounding  quasi-spherical  shells 
of  molecules  occurs  concurrently  to  both  the  nearest  and  next-nearest  neighbor  shells,  but 
with  more  energy  per  molecule  transferred  more  rapidly  to  the  first  shell. 

On  Statistical  Mechanics  of  Small  Systems:  Accurate  Analytical  Equation  of  State  for 
Confined  Eluids,  Mehrdad  Khanpour,  Luis  A.  Rivera-Rivera,  Thomas  D.  Sewell;  Physics 
and  Chemistry  of  Liquids:  An  International  Journal,  Vol.  53  No.  4,  467-480,  2015.  Doi: 
10.1080/00319104.2015.1006631. 

Abstract:  Some  relations  are  derived  using  statistical  mechanics  to  describe  the  effects  of 
surroundings  on  the  properties  of  systems  for  sizes  below  the  thermodynamic  limit.  A 
general  expression  for  the  free  energy  of  closed,  small  systems  is  derived  and  then  used 
to  obtain  the  dependence  of  the  thermal  properties  on  density  and  temperature,  including 
general  expressions  for  equations  of  state  and  internal  energies.  Comparisons  between 
predictions  of  the  current  theory  and  the  results  of  molecular  dynamics  (MD)  simulations 
are  made  for  3D  hard-sphere  and  Lennard-Jones  fluids  for  which  the  surroundings  are 


modelled  as  reflecting  hard  walls  that  confine  the  system  along  one  direction.  The 
analytical  predictions  are  in  excellent  agreement  with  the  MD  results. 

Molecular  dynamics  simulations  of  shock  wave  propagation  across  the  crystal-melt  interface 
of  nitromethane,  Shan  Jiang,  Thomas  D.  Sewell,  and  Donald  L.  Thompson;  manuscript  in 
preparation. 

Abstract:  Molecular  dynamics  (MD)  simulations  were  used  to  study  shock  wave  passage 
with  normal  incidence  across  the  interface  between  (lOO)-oriented  nitromethane  in 
equilibrium  with  the  melt.  The  simulations  were  performed  using  the  fully  flexible,  non¬ 
reactive  SRT  force  field  [Sorescu,  Rice,  and  Thompson,  J.  Phys.  Chem.  B  104,  8406  (2000)]. 
A  reverse-ballistic  geometry  was  used  wherein  the  sample  containing  the  interface  was 
impacted  onto  a  rigid,  stationary  piston  for  the  cases  of  shock  propagation  in  both  directions, 
that  is,  for  shocks  initiated  in  the  crystal  and  transmitted  into  the  liquid  and  vice  versa.  In  all 
cases  the  piston  was  composed  of  the  same  phase  as  the  material  that  impacted  it,  for  example, 
liquid  for  the  case  of  shock  propagation  from  the  liquid  into  the  crystal.  The  local  kinetic 
energies  (intermolecular,  intramolecular,  and  total)  and  stress  states  differ  significantly 
between  the  liquid  and  crystal  regions  and  also  depend  on  whether  the  shock  was  initiated  in 
the  crystal  or  liquid.  The  overall  number  and  spatial  distribution  of  shock-induced  molecular 
orientational  defects  in  the  crystal  for  shocks  initiated  in  the  crystal  are  similar  to  these 
obtained  for  analogous  simulations  for  a  pure  crystal  sample  (no  interface).  By  contrast, 
substantial  differences  in  the  extent  and  distribution  of  shock-induced  orientational  defect 
structures  in  the  crystal  region  were  observed  when  the  shock  was  initiated  in  the  liquid. 
Analysis  of  the  mechanical  states  of  the  system  for  the  two  shock  directions  show  that  the 
damaged  crystal  resulting  from  initiation  of  the  shock  in  the  crystal  is  less  able  to  support 
shear  stresses  than  when  the  shock  is  initiated  in  the  liquid. 


